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Abstract 

A unified theory is advanced to describe both the lateral Goos-Hanchen (GH) effect and the 
transverse Imbert-Fedorov (IF) effect, through representing the vector angular spectrum of a 3- 
dimensional light beam in terms of a 2-form angular spectrum consisting of its 2 orthogonal polar- 
ized components. From this theory, the quantization characteristics of the GH and IF displacements 
are obtained, and the Artmann formula for the GH displacement is derived. It is found that the 
eigenstates of the GH displacement are the 2 orthogonal linear polarizations in this 2-form repre- 
sentation, and the eigenstates of the IF displacement are the 2 orthogonal circular polarizations. 
The theoretical predictions are found to be in agreement with recent experimental results. 

PACS numbers: 41.20.Jb, 42.25.Gy, 42.25.Ja 
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I. INTRODUCTION 

In 1947, Goos and Hanchen [l| experimentally demonstrated that a totally reflected 
light beam at a plane dielectric interface is laterally displaced in the incidence plane from 
the position predicted by geometrical reflection. Artmann Q in the next year advanced a 
formula for this displacement on the basis of a stationary-phase argument. This phenomenon 
is now referred to as Goos-Hanchen (GH) effect. In 1955, Fedorov [3j expected a transverse 
displacement of a totally reflected beam from the fact that an elliptical polarization of the 
incident beam entails a non-vanishing transverse energy flux inside the evanescent wave, 
mbert J4] calculated this displacement using an energy flux argument developed by Renard 
5I for the GH effect and experimentally measured it. This phenomenon is usually called 
Imbert-Fedorov (IF) effect. The investigation of the GH effect has been extended to the cases 
of partial reflection and transmission in transmitting configurations and to other areas 



of physics, such as acoustics [8j, nonlinear optics 



9|, plasma physics [lOj, and quantum 



mechanics (5), And the IF effect has been connected with the angular momentum 

conservation and the Hall effect of light [12I, [ljj. But the comment of Beauregard and 
Imbert [ijj] is still valid up to now that there are, strictly speaking, no completely rigorous 
calculations of the GH or IF displacement. Though the argument of stationary phase was 
used to explain [2| the GH displacement and to calculate the IF displacement 15], it was on 
the basis of the formal properties of the Poynting vector inside the evanescent wave 14] that 
the quantization characteristics were acquired for both the GH and IF displacements in total 
reflection. On the other hand, it has been found that the GH displacement in transmitting 
configurations has nothing to do with the evanescent wave {7]. 

The purpose of this paper is to advance a unified theory for the GH and IF effects 
through representing the vector angular spectrum of a 3-dimensional (3D) light beam in 
terms of a 2-form angular spectrum, consisting of its 2 orthogonal polarizations. From 
this theory, the quantization characteristics of the GH and IF displacements are obtained, 
and the Artmann formula |2j for the GH displacement is derived. The amplitude of the 
2-form angular spectrum describes the polarization state of a beam in such a way that 
the eigenstates of the GH displacement are the 2 orthogonal linear polarizations and the 
eigenstates of the IF displacement are the 2 orthogonal circular polarizations. 
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II. GENERAL THEORY 



Consider a monochromatic 3D light beam in a homogenous and isotropic medium of 
refractive index n that intersects the plane x = 0. In order to have a beam representation 
that can describe the propagation parallel to the x-axis, the vector electric field of the beam 
is expressed in terms of its vector angular spectrum as follows 1g| . 

/*oo roo 

E(r) = — / / A{k y , k g ) exp(zk • r)dk y dk z , (1) 

^ J —oo J —oo 

where time dependence exp(— iujt) is assumed and suppressed, A = ( A x A y A z ) T is the 
vector amplitude of the angular spectrum, k = ( k x k y k z ) T is the wave vector satisfying 
k^ + ky + k^ = k 2 , k = 2mr/X , X is the vacuum wavelength, superscript T means transpose, 
and the beam is supposed to be well collimated so that its angular distribution function is 
sharply peaked around the principal axis (k y o, k z p) and that the integration limits have been 



extended to ±00 for both variables k y and k z 



17| ■ When this beam intersects the plane 



x = 0, the electric field distribution on this plane is thus 

*(y,z) = E(r)U =0 = i- / / Ae^ +k ^dk y dk z , 



hereafter the integration limits will be omitted as such. The position coordinates of the 
centroid of the beam (CQ) on the plane x = are defined by 

= ff&y*dydz = JJ&Hl±dk y dk z 
{V) ff&Vdydz JfA^Adk y dk z [ ' 

and 

ff&zVdydz f f AU&dkydk, 
{Z) JfWmydz JjA^Adk y dk z ' 1 } 

where means partial derivative with respect to k y with k z fixed, means partial 
derivative with respect to k z with k y fixed, and superscript f stands for transpose conjugate. 

Since the Fresnel formula for the amplitude reflection coefficient at a dielectric interface 
depends on whether the incident plane wave is in s or p polarization, it is profitable to 
represent the vector amplitude of the angular spectrum in terms of its s and p polarized 
components. To this end, let us first consider one plane-wave element of the angular spec- 
trum whose wave vector is k° = ( kcos9 ksm8 ) T , where 6 is its incidence angle. Its 
vector amplitude is given by A = A° + A° = A s s° + A p p°, where A s and A p are the 
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complex amplitudes of A° and A°, respectively, s° = ( 1 ) T is the unit vector of A° and 
is perpendicular to the plane xoy, and p° = ( — sin 9 cos 9 ) T is the unit vector of A° and 
is parallel to the plane xoy. This means that A can be represented as 



A 



( -sin0^ 

cos# 

1 



A. 



where 



.4 



A,. 



= A s s + A p p 



(4) 



is what we introduce in this paper and is referred to as the 2-form amplitude of the angular 

(A (° 

spectrum, s = represents the normalized state of s polarization, and p = 

W v. 

represents the normalized state of p polarization, s and p form the orthogonal complete set 
of linear polarizations. 




FIG. 1: Schematic diagram for the rotation of the vector amplitude A of an arbitrarily polarized 
light beam. 



After this element is rotated by angle p around the x-axis as is displayed in Fig. [T], its 
wave vector becomes 



k = M(^)k° 
and its vector amplitude becomes 



^ k cos 9 

k sin 9 cos tp 
y k sin 9 sin p j 



A = M{p)A° = A s s + A p p, 



(5) 
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where 



/ 



M(<p) 



1 

cos <p — sin <f 
y sin cp cos cp j 



is the rotation matrix, 



s = M(»s 

is the unit vector of A s = M(tp)A®, and 



— sin (p cos <p 



( 



\ 



sin i 

p = M(y?)p° = cos 9 COS 

^ cos 9 sin y 

is the unit vector of A p = M(y>)A°. This shows that the vector amplitude (jSJ) can be 
represented as 



{ Px \ 



A 



\ 



Py 

Pz 



A = PA, 



(6) 



/ 



where matrix P represents the projection of 2-form amplitude A onto vector amplitude A 
and is thus referred to as projection matrix, and 

k 7 . 



s y = — sin ip = 

s z = cos <f 

Px = — sin 9 = 

p y = cos 9 cos <p = 

p z = cos 9 sin ip = 



(k 2 y + k 2 z ) 1/2 
k 

h x hy 



k(k 2 y + k 2 z y/ 2 

Now we have successfully represented, through the projection matrix, the vector ampli- 
tude A in terms of the 2-form amplitude A, that is to say, in terms of the 2 orthogonal 
linear polarizations A s and A p . It should be pointed out that in this representation, s is not 
necessarily perpendicular to the plane xoy, and p is not necessarily parallel to this plane. 
Denoting k. r = k r e r = k y e y + k z e z , where k y = k r cosy?, k z = k r simp, e y and e z are the unit 
vectors in y and z directions, respectively, and e r is the unit vector in the radial direction, 



we find that s is in fact the unit vector in the azimuthal direction, s = b v . Furthermore, 
letting p r = p y e y + p z e z , it is apparent that p r = %e r . In other words, p r is in the radial 
direction. The directions of s and p r are schematically shown in Fig. [2j 




FIG. 2: s and p r are in the azimuthal and radial directions, respectively. 

Unit vectors s and p and the wave vector k are orthogonal to each other and thus satisfy 
the following relations, 

4 + s 2 z = i 



Vl + Pi + Vl = 1 



SyPy + S z p z = 

J*hySy | Aj 2 S g 

k x Px + kyPy + k z p z = 



(?) 



The first 3 equations guarantee 



A f A = A ] A. (8) 
From expression ([6]) for the vector amplitude and with the help of Eq. ([7]), we obtain 

(9) 
(10) 



dh y k(k y ~\~ k 2 ^j 



k k 

x — (A* A - A* A ) 



A T — = A T — + 



kxky -(a:a p -a;a s ). 



dk z dk z k(k 2 y + k 2 z ) 

Eqs. (121), ([3]), (E]), (EJ), and ( TTOl) are the central results of this paper, from which the GH and 



IF displacements are derived below. 
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III. DESCRIPTION OF INCIDENT AND REFLECTED BEAMS 



Without loss of generality, we consider an arbitrarily polarized incident beam of the 
following 2-form amplitude, 

A = (Ins + lap) A = LiA, (11) 

where = Ins + l^p describes the polarization state of the beam and is assumed to satisfy 
the normalization condition 

L\~Li = 1, (12) 

angular distribution function A(k y , k z ) is assumed to be a positively-definite sharply-peaked 
symmetric function around the principal axis (kyo,k z o) = (ksin8o,0) and satisfy the nor- 
malization condition 

J J A 2 (k y ,k z )dk y dk z = l, (13) 

and #o stands for the incidence angle of the beam. Eqs. (|12p and ([TBI guarantee the following 
normalization condition for the 2-form amplitude (TTTj) . 




A\Aidk y dk z = 1. (14) 

One example of such a distribution function that satisfies normalization condition ({TBI is 
the following Gaussian function [if], Q] , 



7T 



vol 



exp ( -fkl ) , (15) 



where w y = w / cos9 , w z = w , w is half the width of the beam at waist. A9 = is half 
the divergence angle of the beam. 

According to Eq. ([6]), the vector amplitude of the incident beam is given by Aj = PAj. For 
a uniformly polarized beam that was obtained from a linearly polarized beam in experiments 
19|, 120, |2lJ] , the s components of all its plane- wave elements are in the same direction, and the 
same to the p components. But in our representation advanced here, the s polarizations of 
different plane-wave elements are generally in different directions; so are the p polarizations. 
Considering Eqs. (Q, (iTTj) and (|T5|) together, one concludes that in order to describe a 
uniformly polarized beam mentioned above, it is essential that the incidence angle #o be 
much larger than A6. So we will only consider the case of large 8q below. Fortunately, this 
is just what we have in the case of total reflection. 
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It will be convenient to express Lj on the orthogonal complete set of circular polarizations 
as follows, 

Li = cuf + c i2 i = UCi, (16) 
where en represents the complex amplitude of right circular polarization, q 2 represents 
the complex amplitude of left circular polarization, r = Us = —m\ is the normalized 

\-') 

(A 

state of right circular polarization, I = Up = 4= I is the normalized state of left circular 

w 

1 ( 1 l \ ~ ( c *l 

polarization, U is the unitary transformation matrix U = I I , and Ci = I 

\-i i J \c i2 

f and I form the orthogonal complete set of circular polarizations. Unitary transformation 
guarantees ZJZj = C\Ci. 

When the beam is reflected at plane x — 0, the reflected beam has the following 2- form 
amplitude, 

A r = RiAi = L r A, (17) 

where 

/ R, 



R p 



is tiu > ifiicclioii coefficient matrix. L r R/Lj = | | describes the polaiization stat<^ 

of the reflected beam, and R s = \R S \ exp(i$ s ) and i? p = exp(i$ p ) are the reflection 
coefficients for s and p polarizations, respectively. It will be convenient to express L r on the 
orthogonal complete set of circular polarizations as follows, 

L r = c rl f + c r2 / = UC r , (18) 

where c ri represents the complex amplitude of right circular polarization for reflected beam, 
c r 2 represents the complex amplitude of left circular polarization, C r — I I = R c Cj, 



C r 2 



and R c = U^R/U. Unitary transformation guarantees L\L r = Cj.C r . 



8 



IV. GH EFFECT AND ITS QUANTIZATION 



Applying Eqs. (J2J), ©, and to Ai produces the y coordinate of the centroid of the 
incident beam on the plane x = 0, 

(V)i = o. 

Since R s and R p are all even functions of k z , we have for the y coordinate of the centroid of 
the reflected beam on the plane x = 0, on applying Eqs. (J2]), (JHJ), and ([9]) to A r , 

(V)r --iff (\lrl\ 2 ^ + M^) A*dk y dk z , (19) 



ft J J \ dk y dk y y 

where 

ft = J j \\l rl \ 2 + \l r2 \ 2 )A 2 dk y dk z (20) 
describes the reflectivity of a 3D beam. The above equation can also be written as 

ft = |iti| 2 5H a + l^l 2 ^, 

with 




and 



ft, = I I \R s \'A z dk y dk z 



m p = I I \R p \ 2 A 2 dk y dk z . 




The displacement of (y) r from (y)i is the GH effect and is thus given by 

Don = -^f f \R s \ 2 A 2 ^dk y dk z - M / / \R p \ 2 A 2 ^dk y dk z . (21) 



X J J ' ' Oky ' ^ ft j j ' P ' afcj, 

It is obviously quantized with eigenstates the s and p polarization states. The eigenvalues 
are 

D GHj = ~^j j iRjl'A^dkydh, 

with j = s,p. When the angular distribution function A(k y , k z ) is so sharp that |^ and 
are approximately constant in the area in which A is appreciable, we arrive at the Artmann 
formula Q], 

Dc m = -gt. (22) 

It is now clear that the quantization description of GH displacement depends closely on the 
2-form representation of the angular spectrum. 
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A. Total reflection 



When the beam is totally reflected, the reflection coefficients take the form of 



R s = exp(z$. 



R p = exp(i$ p ), 



(23) 



and 9t = 1. Substituting Eq. fl23|) into Eq. ff2T]) . we obtain 



D 



GH 




hi 



>d$ f 



dk, 



If ff* and -Kjf- are approximately constant in the area in which A is appreciable, the reflected 



beam maintains the shape of the incident beam 22] and the GH displacement takes the form 
of 



D 



GH 



\hi\ 



dk,, ' 



which leads naturally to the Artmann formula (1221) for s or p polarization and agrees well 



with the recent experimental results [19|, [2 



B. Partial reflection and generalized GH displacement 

When the beam is partially reflected, the reflected beam is also displaced from (y)i to 
(y) r in ^/-direction. This is the so-called generalized GH displacement [7| and is given by Eq . 



( I2ip . Such generalized GH displacements may also occur in attenuated total reflection 10(, 



amplified total reflection [23j, and in reflections from absorptive [24( and active [25] media. 



If ||^ and -or- are approximately constant in the area in which A(k y , k z ) is appreciable, Eq. 



211) reduces to 

\lii\ 2 Msd$, |Z i2 | 2 JHp^ 

i->GH = ~ 



m dk y 9t dk y 



which also leads to the Artmann formula (|22l) for s or p polarized beams. 



V. IF EFFECT AND ITS QUANTIZATION 

Now let us pay our attention to the problem of the IF effect. As before, we first want 
to find out the z coordinate of the centroid of the incident beam on the plane x — 0. On 
applying Eqs. (j3J), (JED, and ( jTUl) to Aj and with the help of Eq. (TTBT) . we have 

(^ = (|Q!| 2 -|q 2 | 2 )(^, (24) 
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where 



Wf = Iff ^fv/ ik * ik ' 



Eq. ( 1241) shows that does not vanish and is quantized with eigenstates the 2 circular 
polarizations. The eigenvalues are the same in magnitude and opposite in direction. For the 



Gaussian distribution function (j!5p . we have 26] 



tan 6*o 

at large incidence angle, 9 Q ^> A6 1 . 

The non-vanishing transverse displacement of the incident beam from the plane z = is 
in fact an evidence of the so-called translational inertial spin effect of light that was once 



discussed by Beauregard [23]. Beauregard found that although the transverse wave vector 
of a 2-dimensional beam is identically zero, the 2 circular polarizations have non-vanishing 
transverse Poynting vector, and called this phenomenon the translational inertial spin effect. 
The problem is that the electromagnetic field of so defined 2-dimensional beam is uniform 
in transverse direction, = 0. In order to observe this effect, it is necessary to have a 
bound beam that is not transversely uniform, provided that the expectation of transverse 
wave vector is zero. The 3D beam that we consider here is such a beam satisfying 



(k z ) = A^k z Adk y dk z = A ] k z Adk y dk z = 






FIG. 3: Normalized intensity distributions of right circularly polarized beam (solid curve) and left 
circularly polarized beam (dashed curve) on the z-axis, where O = 10°, A9 = 10" 3 , and z is m 
units of A. 
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For example, when 8q = 10°, we have ~ 0.9A. This displacement has been confirmed 
by the numerical calculation of the field intensity distribution, |\I/(0, z)| 2 , on z-axis as is 
shown in Fig. [31 where the Gaussian distribution function (1151) is considered with A6 = 
1CT 3 <C 9q. The right circularly polarized beam (solid curve) is displaced 0.9A to the 
positive direction, and the left circularly polarized beam (dashed curve) is displaced 0.9A to 
the negative direction. 

When the beam is totally reflected, the 2-form amplitude of the reflected beam is repre- 
sented by Eq. (pH), with R s and R p given by Eq. ([23]). Applying Eqs. ([31. (151). and (flO]) 
to this amplitude and with the help of Eq. ( fTBD gives the transverse displacement of the 

nnnn 

reflected beam from the plane z = 0. So defined displacement is the IF effect |4J, [14J, [19|, [20( 
and is given by 

D IF = (z) r = \J J (M 2 - \c r2 \ 2 )^^A 2 dk y dk z . (25) 

This shows that the IF displacement of the reflected beam is quantized with eigenstates the 
2 circular polarizations. The eigenvalues are the same in magnitude and are opposite in 
direction. Eqs. (1241) and (j25p indicate that the quantization description of IF displacement 
depends closely on the 2-form representation of the angular spectrum. 



In order to compare with the recent experimental results 19|, [20|, [2jJ , we consider such 
an incident beam that has the following elliptical polarization and Gaussian distribution 
function, 

(cosih \ 
A G , (26) 
e sin ip I 

where < ip < n/2. In this case, the IF displacement of totally reflected beam is 

D„ = }VEl^m J J J^ e -*-W e -« smW + 4s - %)dkydK . (27) 

Since Eq. ([27]) holds whether the beam is totally reflected by a single dielectric interface 
[19I or by a thin dielectric film in a resonance configuration [2oj], it is no wonder that the 



observed IF displacement in the resonance configuration 20] is not enhanced in the way 
that the lateral GH displacement is enhanced. 

If the total reflection takes place at a single dielectric interface and the incidence angle 
is far away from the critical angle for total reflection and the angle of grazing incidence 



in comparison with A9, the first and the last factors of the integrand in Eq. ([271) can be 

12 



regarded as constants for a well-collimated beam 22j and thus can be taken out of the 

integral with k y , k z , $ s , and <3> p evaluated at k y = k yQ and k z = k z0 , producing 

sin(2^) sin(0 + $ s0 - % ) (OQ , 
k tan #o 

This shows that for given 8q, the magnitude of Djp is maximum for circularly polarized 
reflected beams (tp = it /4 and <ft + $ s o — &po = (m + 1/2)tc). It also shows that the non- 
vanishing IF displacement for the case of oblique linear polarization of the incident beam 
(0 = m7r) {4] results from the different phase shifts between s and p polarizations in total 
reflection. The incidence angle dependence ~ 1/ tan# is in consistency with the recent 



experimental result [21] . Since #0 is larger than the critica 
wonder that the IF displacement is of the order of Xq/2ti 



angle for total reflection, it is no 



19|,|20|. 



VI. CONCLUDING REMARKS 

We have advanced a unified theory for the GH and IF effects by representing the vector 
angular spectrum of a 3D light beam in terms of a 2-form angular spectrum consisting of 
the s and p polarized components. The 2-form amplitude of the angular spectrum describes 
the polarization state of a beam in such a way that the GH displacement is quantized with 
eigenstates the 2 orthogonal linear polarizations and the IF displacement is quantized with 
eigenstates the 2 orthogonal circular polarizations. We have also derived the Artmann for- 
mula for the GH displacement and found an observable evidence of the so-called translational 
inertial spin effect that was discussed more than 40 years ago It was shown that the IF 
displacement is in fact the translational inertial spin effect happening to the totally reflected 
beam. 

In the 2-form representation of a bound beam presented here, only large incidence angle #0 
in the angular distribution function A(k y , k z ) corresponds to the uniformly polarized beams 
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21| . When 6q is very small, especially when 8q = 0, this representation gives quite 



different beams with peculiar polarization distributions, which needs further investigations. 
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